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Abstract 

We argue that soft non-factorisable terms give a significant contribution to the baryon form 
factors at intermediate momentum transfers and set up a framework for the calculation of 
such terms in the light-cone sum rule approach. Among them, contributions of three-quark 
states with different helicity structure compared to the leading twist prove to be the most 
important. The leading-order sum rules are derived and confronted with the experimental 
data. 
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1 Introduction 


Electromagnetic form factors of the nucleon present a classical observable that characterizes 
the nucleon’s spatial charge and current distributions. The earliest investigations of the 
proton form factor [|I| established the dominance of the one-photon exchange process in the 
electron-proton scattering. The matrix element of the vector current taken between nucleon 
states is conventionally written in terms of the Dirac and Pauli form factors Fi{Q‘^) and 
F 2 {Q‘^)i respectively: 


{P-q\f^{Q)\P) = N{P-q) 




^ 2M 


F2{Q^) 


N{P), 


( 1 . 1 ) 


where P^ is the nucleon four-momentum in the initial state, M is the nucleon mass, P^ = 
(F — qy = M^, is the (outgoing) photon momentum, = —g^, and N{P) is the nucleon 
spinor. The values of Dirac and Pauli form factors at = 0 dehne the electric charge and 
the anomalous magnetic moment of the nucleon: 


Ff (0) = 1, F[^(0) = 0 , F|(0) = Kp = 1.79 , F2"(0) = = -1.91. (1.2) 

Hereafter, ‘p’ and ‘n’ stand for the proton and the neutron, respectively. From the experi¬ 
mental point of view it is more convenient to work with the electric G_e(Q^) and magnetic 
Gm(Q^) Sachs form factors dehned as 

Gm(Q") = Fi(Q=) + F2(Q") , Ge(Q") = F,(Q") + , 

0^,(0)=/*!, = 2.79, G£(0) = =-1.91. (1.3) 

In a special frame of reference, the Breit-frame, Ge(Q^) corresponds to the distribution of 
electric charge and Gm(Q^) to the magnetic current distribution. In the same frame Gm(Q^) 
stands for the helicity conserving amplitude, while Ge(Q^) corresponds to a helicity-flip. 
In the inhnite momentum frame Fi(Q^) and F 2 {Q‘^) are helicity conserving and helicity 
violating, respectively. 

It is known that the experimental data for G'm(Q^) at values of up to 5 GeV^ are 
very well described by the famous dipole formula: 

Tg^(Q2) ~ ~ ( = Gd(Q2) ; i^l ~ 0.71 GeV^ . (1.4) 

fJjp Pn G h/ /hoJ 

For the electric form factor a dipole behavior is observed for below 1 GeV^. For larger 
momentum transfers the experimental situation was unclear until recently, because of SLAG 
0] data contradicting the older DESY results [Q . Both of these measurements were based on 
the traditional Rosenbluth separation of the cross section. Very recently, the Jefferson Lab 
Hall A Gollaboration determined the ratio Gg(Q^)/G^(Q^) from a simultaneous measure¬ 
ment of longitudinal and perpendicular polarization components of the recoil nucleon [^, . 

A systematic deviation from the dipole behavior for the electric form factor was observed 
conhrming the tendency seen in earlier measurements at DESY. 


1 







The theoretical calculation of the form factors from the underlying field theory presents 
a classical problem of the physics of strong interactions. For very large the hard gluon 
exchange contribution proves to be dominant. The corresponding formalism was developed 
in IP, presenting one of the highlights of perturbative QCD. This approach introduces the 
concept of hadron distribution amplitudes as fundamental nonperturbative functions de¬ 
scribing the hadron structure in rare parton conhgurations with a minimum number of Fock 
constituents (and at small transverse separations). The Dirac form factor Fi(Q^) exhibits 
the leading asymptotic 1/Q^ scaling behaviour and can be written in a factorized form as 
a convolution of two nucleon distribution amplitudes of leading twist and a calculable hard 
part. On the other hand, the Pauli form factor F 2 {Q‘^) turns out to be additionally sup¬ 
pressed by an extra power of 1/Q^. It is, therefore, of higher twist and is not accessible 
within the standard approach. 

In practice the calculations of nucleon form factors for realistic values of using the 
standard hard-scattering picture appear to be not convincing. For the asymptotic form of 
the leading twist three-quark distribution amplitude the proton form factor turns out to 
be zero to leading order |^, while the neutron form factor compared to the data is small 
and of opposite sign. A remedy that has been suggested in [§] is that nucleon distribution 
amplitudes at intermediate momentum transfers may deviate strongly from their asymptotic 
shape. Chernyak and Zhitnitsky have shown that it is indeed possible to get a satisfactory 
description of the magnetic form factors using very asymmetric distribution amplitudes in 
which a large fraction of the nucleon momentum is carried by one valence quark. One 
drawback of this suggestion is, however, that even at large values of ~ 10 GeV^ the “hard” 
amplitude is dominated by small gluon virtualities 0, casting doubt on the consistency of 
the perturbative approach. Attempts have been made to increase the region of applicability 
of perturbative QCD by resumming Sudakov-type logarithmic corrections to all orders [|^. 
Unfortunately, the Sudakov suppression of large transverse separations is most likely not 
strong enough to suppress nonperturbative effects, see e.g. 


11 for a detailed discussion. 


Another point of view that is becoming increasingly popular in recent years is that the 
onset of the perturbative QCD regime in exclusive reactions is postponed until very large 
momentum transfers and nonperturbative so-called “soft” or “end-point” contributions to 
exclusive reactions play a dominant role at present energies. In particular, it proves to be 
possible to get a good description of the existing data by the “soft” contribution alone, 
modeled by an overlap of nonperturbative wave functions, see [^]. A weak point of this 
approach is a possibility of double counting, with hard rescattering contributions “hidden” 
in mo del-dependent hadron wave functions. 

In this paper we develop an approach to the calculation of baryon form factors based on 
light-cone sum rules (LCSR) Although the LCSR predictions do involve a certain model 
dependence and the leading-order sum rules may be not very accurate, this technique offers 
an important advantage of being fully consistent with QCD perturbation theory. LCSRs 
reveal that the distinction between “hard” and “soft” contributions appears to be scale- and 
scheme-dependent [Q. It was demonstrated for the case of the pion that the contribution 
of hard rescattering is correctly reproduced in the LCSR approach as a part of the 0{as) 
correction. In recent years there have been numerous applications of LCSRs to mesons, see 
ISl for a review. Baryon form factors, however, were never considered. One reason for this 
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is that the LCSR calculations require certain knowledge about the distribution amplitudes 
of higher twists and for baryons these were not available until recently |^. Another reason 
is that, as we will see, the LCSR formalism for baryons appears to be considerably more 
cumbersome. 

Apart from resolving several technical issues, our main hnding in this work is that the soft 
contribution to the nucleon form factors is dominated by valence quark conhgurations with 
different helicity structure compared to the leading-twist amplitude. Large contributions 
of the distribution amplitudes with “wrong” helicity are important for the electric form 
factor and in a more general context can explain why helicity selection rules in perturbative 
QCD appear to be badly broken in hard exclusive processes at present energies. The sum 
rules in the present paper are derived to leading order in the QCD coupling. Comparing the 
results with the available data we conclude that nucleon distribution amplitudes that deviate 
signihcantly from their asymptotic shape are disfavored. With the asymptotic distribution 
amplitudes, the accuracy of the sum rules proves to be of order 50% in the range 1 < < 

10 GeV^ both for the proton and the neutron. We believe that the accuracy can be improved 
signihcantly by the calculation of 0{as) corrections to the sum rules and especially if lattice 
data on the moments of higher-twist distribution amplitudes become available. 

The presentation is organized as follows. In Section 2 we introduce the necessary notation 
and explain basic ideas and techniques of the LCSR approach on the example of the leading- 
twist contribution to the sum rule. Section 3 contains the derivation of sum rules including 
higher twist corrections, which is our main result. The numerical analysis of the LCSRs 
is carried out in Section 4, together with a summary and discussion. The paper has two 
Appendices devoted to technical aspects of the calculation: In Appendix A we collect the 
necessary expressions for the conformal expansions of nucleon distribution amplitudes and 
in Appendix B a derivation of nucleon mass corrections to the sum rules is given. 

2 Getting Started: Leading Twist 


The method of light-cone sum rules combines the standard technique of QCD sum rules 
r?! with the specihc light-cone kinematics of hard exclusive processes. As the most apparent 


distinction, the short-distance Wilson operator product expansion in contributions of vacuum 
condensates of increasing dimension is replaced by the light-cone expansion in terms of 
distribution amplitudes of increasing twist. Our calculation is similar to the calculation of 
the pion form factor in Ref. [|H . 


Throughout this work we consider the following correlation function 


TAP,q) = i / 


( 2 . 1 ) 


which includes the electromagnetic current 

= euU'juU -F Cdd'^pd 

and an interpolating nucleon (proton) held [§ 

hcz(O) = A^(0)] 75 , 


( 2 . 2 ) 
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Figure 1: The tree-level contribution to the correlation function Eq. (|2.1j). 


(OhcziP) = h{P-z)^N{P). 


(2.3) 


Here z is a light-cone vector, = 0, and the coupling /tv determines the normalization of the 
leading twist proton distribution amplitude 0 . This choice is convenient for our purposes as 
it leads to the same hierarchy of the contributions of different twists as in the perturbative 
approach [^], see below. From the dehnition in Eq. (|1.1|) the contribution of the nucleon 
intermediate state in the correlation function Eq. O is readily derived to be 


z^T,{P,q) = 


Jn 


M2 - P'2 


{P'-z)\ 2F,{Q^) {PP z) - ■ z) / 


where 


+ F2m 



+ ... (2.4) 


P' = P-q (2.5) 

and the dots stand for higher resonances and the continuum. In order to simplify the Lorentz 
structure we have contracted the correlation function with z'^ to get rid of contributions ~ Zy 
that give subdominant contributions on the light-cone. 

On the other hand, at large Euclidean momenta P'^ and q^ = —Q^ the correlation 
function can be calculated in perturbation theory. The leading order contribution is obtained 
from the diagram shown in Fig. |l|. A simple calculation yields 

1 r iqx 

ZyT^{P, ^ y -/T (75 


X 


2ed (0| P^X(0)n^^(0)<(a;) |P) + 4e, (0| P^X(0)w'^(a:)<(0) |P) , (2.6) 


where a ,/?,7 are quark spinor indices. In the light-cone limit 0 the remaining three- 

quark operator sandwiched between the proton state and the vacuum can be written in terms 
of the leading twist nucleon distribution amplitudes 0, H, |T^ 

4 (0| E'=>‘uUa,z)ul(a^z)d%a.^z) |P) = V, ( 75 ^), + A, (/=7'.CIe 

+ Ti ('f7zN\ . (2.7) 
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Each distribution amplitude Vi,Ai and Ti can be represented as 


F(atp-z)= / 


( 2 . 8 ) 


The integration runs over the longitudinal momentum fractions Xi,X 2 -,x^ carried by the 
quarks inside the nucleon with 'Yhi ^ and the integration measure is dehned as 


Vx = J dxidx 2 dx 3 6{xi + X 2 + xs — 1). 


(2.9) 


The normalization is hxed by 


I>xVi{xj,X2yXz;) ^ fziy 


( 2 . 10 ) 


cf. Eq. 


With these dehnitions, we hnd for the contribution in Eq. (|2.6|) 


z,T’^ = 


f XT. a;3Vi(a;i) , ^ 


Vx 


X2Vl{Xi) 

{q - X 2 Py 


2{P ■ zf jiN{P) + ... (2.11) 


where the ellipses stand for contributions that are nonleading in the inhnite momentum 
frame kinematics P —>■ oo,q ~ const, z ~ 1/P. Note that in this case only Vi contributes 
and there is no Lorentz structure ~ {P ■ zY ^ ^ that would give rise to the Pauli form factor 

F2. 

The common idea of QCD sum rules is to match the dispersion representation in Eq. (|2.4|) 
with the QCD calculation at certain “not so large” Euclidian values of the momentum 
P'^ flowing through the nucleon interpolation current. To this end we can rewrite the 
perturbative result in Eq. (p.ll|) in the form of a dispersion relation with a certain spectral 
density 


so 


VT„(P,«)= / ds^k^2(P.j)VA'(i’) + ... . 


( 2 . 12 ) 


Restricting the region of integration to the mass region below the Roper resonance, Sq ~ 
(1.5 GeV)^, one eliminates contributions other than the nucleon |^. If terms of 0{M‘^/Q"^) 
are neglected — which is consistent with twist-3 accuracy — such a representation is obtained 
easily by the substitution s = (1 — xY)QY or s = (1 — X 2 )Q‘^/x 2 for the contribution of 
the d-quark and the n-quark, respectively. The upper bound in the dispersion integral then 
translates into a lower bound in the integral over the corresponding momentum fraction: 
X 2 > QY/{sq -|- or xs > Q^/(so -l- Q^). Finally we follow the usual QCD sum rules 
procedure to use a Borel transformation to convert the power suppression of higher mass 
contributions into an exponential suppression 


1 


-{q — xPy x{s — P'^) 


1 

- exp > , 

X ^ 1 Mg 


(2.13) 
















As the result, a new variable — the Borel parameter Mb — enters instead of P'^. Equating 
the Borel transformed versions of Eq. (p.4|) and Eq. ( p.ll| ) we finally arrive at the sum rules 


^tw-d(g.) _ 


fN _ 


Cd I VxVi{xi)exp 
2e„ / VxVi{xi)exp 


X 3 MI 

X2Q‘^ — x^M"^ 


X 2 M 


B 


0 1 ^X 3 - 
Q (x 2 - 


+ So 

+ So 


F*'^-3(Q2) = 0. 


(2.14) 


where the superscript ‘tw-3’ reminds of twist-three accuracy^. 

The ©-functions in Eq. ( p.l4 ) originate from the restriction of the spectral density to the 
duality region s < sq and conhne the integration region to values of X 3 —> 1 and X 2 —> 1 at 

^ 00 for the contributions of d and u quarks, respectively. It follows that in our approx¬ 
imation the form factor is dominated by parton configurations where the scattered quark 
carries almost all the momentum of the nucleon. This is precisely the soft, or Feynman mech¬ 
anism, and the standard wisdom tells that this contribution has to be subleading at very large 
momentum transfers. Indeed, using the asymptotic wave function Vi(xi) = 120 /ArXiX 2 X 3 and 
expanding in 1 /Q^ we hnd that the sum rule result in the limit ^ 00 behaves as 


r-so 


^tw- 3 (gg ^ 20 (e, + 2 e„) 


ds 




(2.15) 


i.e. it is suppressed by two additional powers of l/Q^ compared with the expected asymptotic 
behavior. This observation is in agreement with an analysis of the soft contribution in the 

|I 8 | : The soft overlap contribution of two 
asymptotically vanishes as 1 /Q®. 


framework of the Drell-Yan description of Fi 
leading twist light-cone wave functions computed in |]I9 


We will find, however, that this strong suppression does not hold for contributions of wave 
functions with different helicity structure. 

It is well known that the form factor is strongly sensitive to the shape of the nucleon 
distribution amplitude. The general QCD description of distribution amplitudes is based 
on the conformal expansion B 0 0 . There exist several concrete models that take into 
account the first few conformal partial waves. In Fig. ^jwe show the twist-3 LCSR prediction 
(Eq. ( |2.14|) ) for the magnetic form factor (to this accuracy Gm = Fi) normalized to 
the phenomenological dipole parametrization Eq. (|1.4|) , for three different choices of the 
distribution amplitude 

V^^{Xi,fi ~ IGeV) = 120 XiX 2 X 3 /Ar , 

Vi(xi,/i ~ IGeV) = 120 xiX 2 X 3 /Ar 1 ^3 (/i)(l - X 3 ) 

~ IGeV) = 120 xiX 2 X 3 /Ar 1 + ^3 (/i)(l - X 3 ) 


+ <pf[3 - 21 x 3 + 28x2] + [5(x 2 + x2) - 3(1 - X3)2)] ,(2.16) 


*Throughout this paper we give explicit expressions for the proton. Neutron form factors are obtained 
by the obvious substitution e„ e^. 
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Figure 2: Twist-3 approximation to the soft contribution to the magnetic form factor of the proton 
(left) and the neutron (right), see Eq. (|2.14|) . Plotted is the ratio G^“/(/Up/nG£)) where Gd is the 
dipole fit Eq. ([1.41) . The three curves correspond to different models of the leading twist distribution 
amplitudes: The full line corresponds to the dashed line to Vi, and the dashed-dotted line to 
respectively, see Eq. ( 2.16 ). Sum rule parameters are taken to be = 2GeV^ and sq = (1.5 GeV)^. 
The data points are taken from |20, p]]. 


where 

/^(/x = IGeV) = (5.3 =E 0.5) x IQ-^GeV^ , 

0;|'(/x = IGeV) = 1.1 ±0.3, 

(/X = IGeV) = 0.615, 

(/X = IGeV) = 3.68. (2.17) 


The hrst expression in Eq. ( p.l6| ) dehnes the asymptotic distribution amplitude in the limit 
> cx) and corresponds to the contribution of the lowest conformal spin (“S-wave”) 0]. 
The second and the third expressions correspond to taking into account the next-to-lowest 
(“P-wave”) and also the next-to-next-lowest (“D-wave”) conformal spin, respectivelj|]. The 
parameters in Eq. (|2.17|) are taken from the last reference in [^. They were obtained using 
QGD sum rules for the lowest moments of the distribution amplitudes and their accuracy is 
subject of debate. Note that the dependence on /at in fact cancels out of the sum rule in 
Eq. ( p.l4|) . The numerical impact of the scale dependence of the remaining parameters is 
negligible compared to the intrinsic uncertainty of their estimates, so that we do not take 
this effect into account. 

From Fig. ^ it is seen that the results vary by factor 2-3 depending on the choice of the 
distribution amplitude. The “P-wave” approximation in the second line in Eq. (|2.16|) seems 
to £t best while for the asymptotic distribution amplitude the neutron form factor turns out 
to be zero. One has to keep in mind, however, that in this calculation F 2 {Q‘^) is zero both 
for proton and neutron and for arbitrary choice of the distribution amplitude. Hence we 
have, e.g. Gm{Q^) = Ge{Q^) identically, which is clearly not a satisfactory approximation 
in a few GeV region. We conclude that the leading twist approximation to the LGSRs is not 


^We have rewritten the original expressions given in Q using the conformal basis defined in the work 
In this basis, all terms are mutually orthogonal with respect to the conformal weight function V“^(a:i)- 
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sufficient for the quantitative analysis and in the next section we proceed to the complete 
treatment, including higher-twist effects. 


3 Beyond the Leading Twist 


Apart from radiative corrections, the correlation function in Eq. ( p.l|) receives further con¬ 
tributions of higher-twist distribution amplitudes. These contributions arise in two ways. 
First, one has to remember that the hard quark propagator ~ in Eq. (|2.6|) receives cor¬ 

rections in the background color held which are proportional to the gluon held strength 
tensor and give rise to four-particle (and hve-particle) nucleon distribution amplitudes cor¬ 
responding to parton states with a gluon or quark-antiquark pair in addition to the three 
valence quarks. Such corrections are usually not expected to play any signihcant role (see 
e.g. [^), and in this work we do not take them into account. Second, one has to improve 
on the treatment of the matrix elements of the three-quark operators in Eq. (|2.6|) taking 
into account contributions of other Lorentz structures and less singular contributions on the 
light-cone, beyond the leading-twist approximation in Eq. ( p.7|) . 

We only need to retain vector-like Dirac structures since it is easy to see (cf. Eq. (H) 
that all others do not contribute to the sum rule we are considering. To this accuracy, the 
complete decomposition of the three-quark operator reads 


4 (0| M^(aix)M^(a2x)ci^(a3x) \P) = 

= Vi + 


(3.1) 


K (75iV), + V 2 M (i^75iV), + V 3 M ( 7 ,C)„^ 


+V 4 M 2 (75iV)^ + (7^C)„^ (za'^^x,75iV), + 


The expansion in Eq. (|3.1|) has to be viewed as an operator product expansion to the leading 
order in the strong coupling. Each of the functions Vi depends on the deviation from the 
light-cone at most logarithmically, and we have also included the O(x^) correction to the 
leading-twist-3 structure, denoted by V^. The invariant functions Vi(xi),..., V6(xj) can be 
expressed readily in terms of higher-twist nucleon distribution amplitudes Vi(xi),..., V6(xi) 


introduced in Ref. 16 


Vi = El, 2p- XV 2 = Vi-V2-V3, 

2 V3 = R3 , 4 p - XV4 = -2W + R3 + 14 + 2R5 , 

4 p ■ X Vs = R4 - , ( 2 p ■ X ) 2 Ve = - Vi + V2 + V3 + V4 + Vs - Ve . 


(3.2) 


In difference to the “calligraphic” functions Vi(xi),..., V6(xj) each of the distribution am¬ 
plitudes Vi(xi),..., Velxi) has dehnite twist, see Table |I], and corresponds to the matrix 
element of a (renormalized) three-quark operator with exactly light-like separations x^ —> 0, 
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see Table 2 and Appendix C in Ref. 
plitudes V 2 (xj), • • •, V^i^Xi) correspond to 


for the details. The higher-twist distribution am- 
“wrong” components of the quark spinors and have 


















twist-3 

twist-4 

twist-5 

twist-6 

Ri 

R 2 , R 3 

R 4 , R 5 

Re 


Table 1: Twist classification of the distribution amplitudes in Eq. ( p.2|) . 


different helicity structure compared to the leading twist amplitude. For baryons these 
components cannot all be traded for gluons as in the case of mesons [p4| 


bad” 

They are not all 

independent, but related to each other by the exact QCD equations of motion. As the result, 
to the leading conformal spin accuracy the hve functions V 2 {xi ),..., V^lxi) involve only one 
single nonperturbative higher twist parameter. In the calculations presented below we use 
the conformal expansions of higher twist distribution amplitudes to the next-to-leading or¬ 
der (include “P-wave”). This accuracy is consistent with neglecting multiparton components 
with extra gluons (quark-antiquark pairs) that are of yet higher spin. Explicit expressions 
for the distribution amplitudes are collected in Appendix A. 

The invariant function corresponding to the corrections of order 0(x'^) to the leading 
-twist Lorentz structure in Eq. ( p.lD is twist-5 and in general is a complicated function of 
distribution amplitudes. In the present context there is an important simplihcation that two 
of the quarks in Eq. (p^) always appear to be at the same space-time point, cf. Fig. |l]. In 
this specihc conhguration, can be obtained using the technique of Ref. that has been 


designed for the studies of mesonic operators. The corresponding calculation is carried out 
in Appendix B and the expressions for in the above mentioned limits present one of the 
new results of this paper. Substituting the full expansion Eq. (|0|) in Eq. (|2.6| ) we obtain 
the result: 


z,T'^ = 2{P ■ zf ^N{P){ed 


Vx 


x:iVi{xi) 

{q - x^py 


(3.3) 


-I- M / dx3 Xg 

Jo 


^ -V2-V3-VA-V5 -) ( 0 : 3 ) 


(g - x^py 


2M^ / dx'iX^ 


1 lRi-R 2 -R 3 -V^ 4 -R 5 + n 

3^3V- 




{q - x^Pf 


+ 2e^ 


X 3 ^ X 2 




M 


— / dxg xl 


1 [V2 + V3 -V,){x,) 


R1-R2-R3-R4-R5 + Re) (X 3 )n 

{q - XsPy 


+ 2e^ 


X 3 <->■ X 2 


+ 


We see that at this stage a contribution to F 2 (Q^) arises. The functions and 

(the second one appears in the u-quark contribution) are given in Eq. (|B.15|) and Eq. (|B.16|), 
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while the distribution amplitudes with a ‘tilde’ are dehned as 


V{x,) 


V(X2) 




^{X2) 




<lxiV{xi,l- Xi -x'^,x '^), 



'‘I—Xo 


dxiV[xi, x' 2 ,1 — xi — x' 2 ), 




'‘I—Xq 


dxiV{xi, 1 - xi - xl,xl ), 






dxiV{xi,X2,1 — xi — X 2 ) 


(3.4) 


and result from partial integration in X 2 or 0 : 3 , respectively: The integration by parts in X 2 
or X 3 is done in order to eliminate the 1/P ■ x factors that appear in Eq. 0 when the 
distribution amplitudes Eq. (|3.2|) are inserted. After this, the / d^x integration becomes 
trivial. The surface terms sum up to zero. 

The Borel transformation and the continuum subtraction are performed by using the 
following substitution rules: 


dx 


q{x) 


dx p(x) 


(g - xpy Jo x2 (s - 

1 dx . . 

-^g{x) exp 
x^ 


MI 4 


xQ2 


xMl 


xM 2\ p(xo)e ^o/M| 

Hf) ^ Q^ + xlAP ' 


(3.5) 


dx 


p(x) 


{q — xP)® 


dx g{x) 

/o x^ (s — P'^)3 

1 dx 


2M' 


B J XQ 


Xn 


X'" 


p(x) exp 


d 




2 g2 + xlM^ 


xMl 

g{xo) 


xM^ 

14 


p(xo) e 


-so/M| 


2 xo (Q2 + a;2M2) Ml 


dxo Xo (Q^ + xlM‘^) 


^-solMl 


(3.6) 


where in difference to Eq. ( p.l4|) we have to keep the nucleon mass in the kinematical relation 
s = + (1 ~ x)M 2 and xq is the solution of the corresponding quadratic equation for 

s = So: 


Xo = 


+ So - M2)2 + 4M2g2 - (g2 + So _ M^)]/{2M^) . 


(3.7) 


The contributions ~ e in Eq. (|3.5|) and Eq. ( p.6|) correspond to the “surface terms” 

arising from successive partial integrations to reduce the power in the denominators {q — 
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xPy^ = (s — with n > 1 to the usual dispersion representation with the 

denominator ~ (s — P'^). Without continuum subtraction, i.e. in the limit sq —>■ C)0 these 
terms vanish. 

Collecting everything together, we get the hnal expressions: 


pf(g2) = 


ed 

In 


dxa exp 

m2 




X 3 MI Ml 




dXiVi{Xi,l-Xi-X^,X:i) 


Ml 

M^ 

Ml 


2 V 1 -C 2 -C 3 -C 4 -W - -Vf ) (X 3 ) 

3^3 


M2(a;0)2 
g2 + {xlfM'^'' 
m2 


Vi-V2-V^-V^-V^ + V^) (X 3 ) 

(^21/i-C2-C3-C4-C5 - (^ 3 ) 




+ ^ ( C1-I/2-C3-C4-I/5 + n ) (x°) 


m2 


d UA?[Vi-V 2 -V^-V^-V^ + V,){xI) 


dx^ 


g2 + (a;0)2M2 




and 

p|(g2) = 


ed dx3 


fN J Xg ^3 


exp 


X^Q"^ a;3M2 


:r3M' 


M2 


[2M2 

(~ ~ ~ \ 

[mi 

(E2 + E3-Eij 


2M^ 


2M^xl 


g2 + (x^)2M2 


Cl-C2-C3-1^4-l>5 + n (X|i) 


m2 


Ml 


( 3 . 8 ) 


^2 + ^ 3 -K C1-C2-C3-C4-C5 + C6 (xlj) 


+ M2x°-^ 
■^dx9 




^ /xO C 1 -I/ 2 -C 3 -I/ 4 -C 5 + C 6 (xO) 


g2 + (xO)2M2 


-{so-M^)/Ml 


( 3 . 9 ) 


The sum rules in Eq. (^) and Eq. (| 3 . 9 |) present the main result of this paper. 

In the limit —>■ 00 the sum rules can be simplihed. Employing asymptotic wave 
functions we get the asymptotic behavior 

^ M2 


pp(g2) = 


37 ^Ai 

y+ 2— )ed + 


3 Jn 


Ml 
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F|(Q2) = _2 


1 / m2 


^ 

^'+15 M2 


ed 


, 8 m2 Ai / 2 m2 

^ “ 9 ii^ + V + 15 ii^ 


M2 

M| 


fSO 


X ( siiM"e-‘*»-"’>/"S + j dss=e-‘*-"’l'" 


(3.10) 


Note that the contribution ~ 1/Q® to Fi(Q 2) arises (cf. Eq. (|2.15|) ), which has two sources: 
The constants ~ 37/3 originate from the nucleon mass corrections and the terms ^i/J n 
correspond to the contributions of higher-twist three-quark operators. Numerically Xi/Jn ~ 
—5, see the next Section, so that both contributions appear to be of the same order. To 
avoid misunderstanding note that to leading order in we are dealing with the soft, or 
Feynman contribution to the form factors only. The true asymptotic behavior Fi 1/Q" 
can be reproduced in the LCSR approach as a part of the 0(«2) correction, a calculation 
which goes far beyond the tasks of this work, see Ref. for the detailed discussion for the 
pion. 


4 Numerical Results and Discussion 


For the numerical evaluation of the form factors we used the standard value of the contin¬ 
uum threshold Sq = 2.25 GeV2 and varied the Borel parameter in the range 1.5 < M^ < 
2.5 GeV2. The results appear to be rather stable so that in the Figures given below we 
take M^ = 2 GeV2 as our standard choice. The solid curves marked ‘asy’ in Fig. and 
Fig. I are obtained using the set of asymptotic distribution amplitudes corresponding to the 
contributions of operators with the lowest conformal spin. As already discussed in Sec. 2, 
the dependence on the normalization of the leading-twist amplitude actually drops out of 
the sum rules so that to this accuracy the form factors depend on a single nonperturbative 
parameter - the ratio of matrix elements of the twist-4 and twist-3 operators 


In 


= -5.1 ±1.7. 


(4.1) 


The constant fjy is dehned in Eq. 
so-called Ioffe current |2^: 


while Ai is the familiar nucleon coupling to the 


miO) = [u\0)Cj^u^{0)] j.rdHO) , 

{0\vi\P) = AiM iV(P). 


(4.2) 


iGl and 


The number in Eq. ( [4.1|) is the QGD sum rule estimate at the scale 1 GeV (see 
references therein). Notice that the Ioffe coupling Ai is large compared to f^. Hence the 
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Figure 3: The LCSR prediction for the soft contribution to the magnetic form factor of the proton (left) 
and the neutron (right). Plotted is the ratio where Gd is the dipole fit Eq. (|1.4D . The 

solid curves marked ‘ASV are obtained using the set of asymptotic distribution amplitudes corresponding 
to the contributions of operators with the lowest conformal spin. The curves marked S/?show the results 
including the next-to-lea ding terms in the conformal expansion with parameters estimated using QCD 
sum rules. The dashed curves are obtained using asymptotic distribution amplitudes and a reduced 
relative normalization of the higher-twist contributions within the error range, see text. The data points 
are taken from \2C, ^|. 




Figure 4: The LCSR prediction for the electric form factor. Plotted is the ratio 

proton (left) and {G'^/GdY the neutron (right), where Gd is the dipole formula Eq. (11.4] ). The 

identification of the curves is the same as in Fig. The data points are taken from |Q, |^, |21j]. 


contribution of three-quark states with a different helicity structure compared to leading 
twist is numerically most important for the form factors in a few GeV range. 

Beyond the leading order in the conformal expansion the parameters become more nu¬ 
merous and less known. The curves in Fig. |]and Fig. ^marked ‘SR’ are obtained including 
the contributions of the next-to-the leading spin (“P-wave”) and using QCD sum rule esti¬ 
mates p, 1^ for the four additional parameters that enter the sum rules to this accuracy. 
For the leading twist this approximation corresponds to the expression in the second line 
in Eq. ( p.lbj) , and cf)'^ in Eq. (|2.17|) is the hrst of the additional parameters. Among the 
remaining three parameters one is of twist-3 and two are of twist-4, see Appendix A for the 
details. It is seen that with large “P-wave” contributions the description of the form factors 
becomes much worse compared to the asymptotic distributions. This is not unexpected, 
since from the experience of calculations of the pion form factor it is known that QCD sum 
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rules based on the local operator product expansion tend to overestimate the higher spin 
corrections considerably, see e.g. p6|. On this evidence, we conclude that large corrections 


to the asymptotic distribution amplitudes of higher twists are unlikely. 

The agreement with the experimental data generally becomes better if the ratio Xi/Jn 
is chosen to be 30% lower compared to the estimate in Eq. ([4.11) , which corresponds to 
the bottom of the given error range. The corresponding results are shown in in Fig. ^ and 
Fig. §by the dashed curves. At the same time, we have checked that the agreement cannot 
be improved substantially by adding “P-wave” contributions to the distribution amplitudes 
with their normalizations taken as free parameters. This, again, can be considered as an 
argument against large corrections to the asymptotic distributions. On the other hand, we 
believe that quantitative conclusions can only be made after the calculation of radiative 
0 {as) corrections to the sum rules, which goes beyond the tasks of this paper. 

To summarize, in this work we have set up a framework for the calculation of baryon 
electromagnetic form factors in the light-cone sum rule approach. The light-cone sum rules 
are derived to leading order in the strong coupling and confronted with the experimental 
data. We argue that soft non-factorisable terms give a signihcant contribution to the baryon 
form factors at intermediate momentum transfers. Among them, contributions of three- 
quark states with different helicity structure compared to the leading twist prove to be the 
most important. A quantitative analysis of baryon form factors requires the calculation of 
radiative corrections to the sum rule, which is a challenging but doable task. The approach 
can easily be generalized to the study of transition form factors like A — N'j and to weak 
decays of heavy baryons. 


Appendices 

A Summary of Nucleon Distribution Amplitudes 


To make the paper self-contained we collect here the necessary information on the set of 
nucleon distribution amplitudes that enter the LCSRs in Eq. (^.8| ) and Eq. ( ^.8| ) for the 
nucleon form factors. For dehniteness, we consider the proton distribution amplitudes. The 
presentation in this Appendix follows Ref. 
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The nucleon distribution amplitudes can be classihed according to twist in the light- 
cone quantization approach of |^. Hereby the quark held operators are decomposed in 
“good” and “bad” or “plus” and “minus” components, respectively: q = + q~. To this 

end it is useful to introduce two light-like vectors p = p+,z = Z-, p^ = z"^ = 0, where 
P/i = Pfi + 0{M‘^). The projection operators on the “plus” and “minus” components are 
A+ {2pz) and A_ =/^/(2pz), respectively. The leading twist-3 amplitude is identihed 

as the one containing three “plus” quark helds while each “minus” component introduces an 
additional unit of twist. This classihcation is adapted to the applications in hard reactions: 
Higher twist amplitudes in general produce power-suppressed corrections to physical cross 
sections. 

The distribution amplitudes Ri(xi),..., in Eq. (|3.2|) can be dehned 0, I, 0 as 

nucleon-to-vacuum transition matrix elements of non-local operators built of quark helds 
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with definite helicity 


^TU) ^ 1(1 


(A.l) 


and separated by light-like distances. In particular, the leading twist-3 distribution am¬ 
plitude Vi corresponding to the (u~^u~^d~^) component of the three-quark operator can be 
written as 


Jjk 


u 


i{aiz)C j^u,{a 2 z)\ 0j,{a^z) \P) = --pz ^nUvx Vi(xi). (A.2) 

J |io) / J 


Here and below we use the notation 


Ui{aiz)C^Uj{a2z) = - {ul{aiz)C^u\{a2z) + ul{aiz)C/u]{a2z) 

J |io) 2 V ■' 

There exist two twist-4 distributions: 


(A.3) 


Jjk 


U. 


(aiz)C/u^(a 2 z) i^dlia^z) jP) = --pz Vxe , 



(01 £ 


ijk 


Ui{aiz)C ^''y±^Uj{a 2 z) ^ \P) = —pzMJVx e 


(A.4) 


corresponding to the {u~^u~^d~) and {u~^u~d~^) projections, respectively. Here T stands for 
the projection transverse to 2 ; and p, e.g. = y'^fi'l^y^ with — {p^z^P z^Pi,)/pz. 

Similarly, the two possible projections with two “minus” fields: {u~u~d'^) and {u~u'^d~) 
give rise to two distribution amplitudes of twist-5 


( 0 | P^'^^Ui{aiz)C ^Uj{a2z)\ / 4 (« 3 ^) \P) = jvx 1 / 5 ( 2 ;,), 


u^{aiz)Cfhf±^uJa 2 z) y^^ 4 (® 3 ^) 1-^) = —pzM']^N^ / Vxe *^^^"^*“* 1 / 4 ( 2 ;^), 
L J |io> J 


(A.5) 


respectively. Finally, there exists a single twist-6 three-quark distribution amplitude 


Jjk 


Ui{aiz)C^uJa2z)] 1^) = ^6(a;,) (A.6) 


!io> 


corresponding to the {u~u~d~) projection. 

The distribution amplitudes are scale-dependent and can be expanded in contributions of 
conformal operators that do not mix with each other under the renormalization (to one-loop 
accuracy). To the next-to-leading conformal spin accuracy the expansion reads [|l6l : 

Vi{xi,p) = 12QX1X2X2, [(j^lip) + 0 ;]'(p)(i - 32:3)], 
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V2{xi,fi) = 24 X 1 X 2 [04(/^) + 04 (/^)(1 - Sa^s)] , 

Vsixi,^) = 12X3 [04(/^)(l - ^3) + 04 (/^)(^1 +xl- 0:3(1 - X3)) + 04^(/i)(l-X3-lOa:ia;2)] , 
14 ( 0 :*,/i) = 3 [ 05 (/i)(l-X 3 )+ 05 '(/i) ( 20 : 10 : 2 -X3(l-X3))+0^(/i)(l-X3-2(a:i+ 0 : 2 ))] , 
V^{xi,ii) = Qxs [(plifi) + 0^(/i)(l - 20:3)] , 

Ve{xi,fi) = 2 [0|!(/i) + 0f[(/i)(l-3o:3)] , (A.7) 


so that it involves 14 parameters. Not all of them, however, are independent; They are 
related to each other by exact QCD eqnations of motion. 

in 


The corresponding analysis was carried ont 


m 


It tnrns ont that the conformal 
expansion coefficients in Eq. ( |A.7| ) can all be expressed in terms of 6 independent matrix 
elements of local operators. To the leading conformal spin accuracy (“S-wave”) only two 
parameters enter, and Ai, which were dehned in Eq. (p.3|) and Eq. ([4.2|) , respectively. 
One obtains 116 


^l = <Pl = fN, 


0 — 05—2 


04 = 05 = 2 


Ai) . (A.8) 


To the “P-wave” accuracy four new parameters enter, corresponding to matrix elements of 
the local operators similar to Eq. (p7^ {Vf, A“) and Eq. (/0 /“) but with an additional 
covariant derivative. The parameters Vi, A“ are leading twist-3 and, in particular, the 
coefficient (fi) appearing in the hrst line in Eq. ( A-7|) is given by 

= 4>t/h = (7/2)(l - 311“). (A.9) 


The parameter describes the component in the leading twist distribution amplitude that 
is antisymmetric in xi ^ X 2 and does not contribute to the sum rules directly. It does 
produce, however, a contribution to the symmetric parts of higher twist amplitudes. In 
addition, there are two new parameters ff, /“ that are genuine twist-4. One obtains for 
twist-4: 

= I (Ai(l - - 4/r) + h(2A'; - 1)) , 

^ = 1 (Ai(3 - lO/f) - /A,(10n'' - 3)) , 

V. 4 - = - j (Ai(2 - Iff + /“) + fniAf + 314“ _ 2)) , 

i’t = (Ai(-2 + 5/,“ + 5/r) + /k(2 + 541 - 511“)) , (A.IO) 

for twist-5: 

C = I (Ai(/? - /,") + /n(241 - 1)) , 

■At = (Ai(4/i“ - 1) + /«(3 + 4V1“)) , 
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Ai/ In 

Vf 


ft 

n 

QCDSR 

-5.1 ±1.7 

0.23 ±0.03 

0.38 ±0.15 

0.6 ±0.2 

0.22 ±0.15 

asymptotic 

-5.1 ±1.7 

1/3 

0 

3/10 

1/10 


Table A: Numerical values for the expansion parameters of the nucleon distribution amplitudes. 

= 5(Ai(/?-/r) + /«(2 -A'!- SV'/)) , 

'Pi = -g bi(-l + 2/? + 2/r) + /iv (5 + 2A'; - 2Vi^)) , (A.ll) 

and for twist- 6 : 

4-^ = i (Ai(l - 4/f - 2/,") + /„(!+ iAX)) , 

<l>t = -i (Ai(l - 2/f) + - 1)) , (A,12) 

Note that there is a mismatch between the twist classihcation of the distribntion ampli- 
tndes that implies connting powers of the large momentum and the twist classihcation 
of the local operator matrix elements: E.g. the parameters of the twist-4 distribution am¬ 
plitudes depend both on the the leading twist-3 matrix elements /at, Vf, A“ and the twist-4 
matrix elements /“. This “propagation” of lower-twist matrix elements to higher- 

twist distribution amplitudes is well known and usually referred to as Wandzura-Wilczek 
contribution. Note that the distribution amplitudes of twist-5 and twist -6 are entirely of 
Wandzura-Wilczek type since there exist no geniune twist-5 and twist -6 operators to this 
order in the conformal expansion. QCD sum rule estimates for the twist-3 and twist-4 
parameters are given in Table ^ together with their asymptotic values in the ^ oo 
limit. 



B The 0(x‘^) Corrections 


In this Appendix we present the calculation of the 0{x‘^) corrections to the light-cone expan¬ 
sion of the three-quark operator in Eq. (|3.1| ). We consider the vector Lorentz projection as 
it is the only one that contributes to the sum rules. We also make use of the important sim- 
plihcation that to leading order in the strong coupling only the bi-local operator with either 
d or u quark shifted from the zero space-time point enters, as can be seen from Eq. (|2.6|) . 
Consider Erst the d-quark contribution: 


[u^C^aU^] iO)d^^ix) \P) 


-X 


(Vi + + V2MP„ 


+ VsM (7a75A^)., + V.M^Xo^ 


(B.l) 
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We remind that Vi starts at leading twist-3, and hence is of twist-5. We do not show 

the O(x^) corrections to the other Lorentz structures since they are yet higher twist and will 
be omitted in what follows]]. The invariant functions Vi,..., Ve can easily be related to the 
distribution amplitudes of dehnite twist by taking the light-cone limit of Eq. (B , ^ 0, 
see Eq. (|3.2|) |]T]. On the other hand, the calculation of is not that immediate since 


in the light-cone limit this contribution vanishes. The meaning of the separation of Vi and 
jg easily understood upon the short distance expansion —>• 0. In this way, the 

nonlocal “string” operator in the l.h.s. of Eq. ([B.l|) is Taylor-expanded in a series of local 
operators with three quark helds and the increasing number of (covariant) derivatives acting 
on the d-quark. The separation of the leading twist part of each local operator corresponds 
to the symmetrisation over all Lorentz indices and the subtraction of traces. Without loss 
of generality, we can consider the matrix element contracted with an additional factor Xa, 
see Eq. o. so that the symmetrization is achieved. To subtract the traces, we formally 
write 


OO 


Xad{x\ = 2 ^^ 
n=0 


• • • vO 




X 

T 


n -I- 1 


l-l'i il^j 


{Xa ■ ■ ■ 




(B.2) 

where ‘It’ stands for the leading-twist part. Observing that the subtracted 

contributions O(x^) can be reassembled in the form of a non-local string operator: 

[uX^u^] {o)d’;{x) \p) = (o| [uX^u^] {o)d’;{x)\_^ \p) 


x‘^ 


dt 


dxr,dx°‘ 


(0| [uX (O)d^(te) \P) . (B.3) 


The same result can be obtained by observing that the leading-twist nonlocal operator 
has to satisfy the homogenous Laplace equation 

fp 

(01 [£«‘ [u<C (0)ci?(i)] „ |F> = 0, (B.4) 

Using QCD equations of motion the second line in Eq. (|B.3|) can be simplified to 

fp 

\uX^un mdtitx) = (0)D„d^(te) + gluons 

dXaOX°‘ L J / L J r 

= 2tdae^^^ {d)d^{tx) + gluons, (B.5) 

where da is a derivative with respect to the overall translation for the matrix element 
we can make the substitution —>■ —iPa- Inserting this result in Eq. (|B.3| ) we hnally obtain 

[uX^u^] {d)d^^{x) \P) = (0| [uX^u^] (0)d^(x)]^jP) 


+ ^H 2 P„) 


dtt (0| (O)d^(fa;) |P) -|- gluons . (B.6) 


^Strictly speaking, since we are not taking into account twist-6 contributions induced by 0{x^) corrections 
to V 2 and V 3 , in order to be consistent we have to discard the contribution of Ve altogether. This contribution 
to the sum rules appears to be numerically neglible, however. 
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Notice that the r.h.s. only involves (up to corrections with additional gluons) the already 
known distribution amplitudes. This equation therefore allows us to determine — 

which appears on the l.h.s. of Eq. ([B. 6 |) — up to gluonic corrections. 

First, consider the leading twist contribution. We can write 


(0| [u^C^u^] \P) = - Px Vi (ysiV). 


' Vx 

Vx 


-iP-xxs 




P ■ xVaM 


V 3 M+... (B.7) 


where [e and [e are the leading-twist components for the free 

helds, dehned as the solutions of the corresponding homogeneous Laplace equation. Note 
that the factor P ■ x in the second line in Eq. (|B.7| ) is not included under the [.. .]it bracket 
since (P • x)V 2 = l/2(Vi — V 2 — V 3 ) is a function of momentum fractions only and does not 
contain any dependence on the position vector x. The solution can easily be constructed 
order by order in the expansion. To the required 0{x‘^) accuracy, the result reads 






dt e 


-iP-xxzt 


(B. 8 ) 




^-iP-xx^ 






dt t e 


2 -iP-xx^t 


+ i . 


X3X 


dt te 


-iP-xx^t 


The corresponding contribution to is proportional to the nucleon mass squared and 

involves the leading twist distribution amplitude, being an exact analogue of the Nachtmann 
power suppressed correction in deep inelastic scattering. The second contribution on the 
r.h.s. in Eq. (|B. 6 |) is special for the exclusive kinematics since it involves a derivative over 
the total translation that vanishes for forward matrix elements. Its explicit form is easily 
found by contracting the three-quark matrix element in Eq. ([B.l|) with Pa instead of Xq, and 
inserting the resulting expression in Eq. (|B. 6 |). One gets 


X 


[-i2Pa) / du(0|p^'^ (O)d^(ta;) |P) 

Jo 


i Vx + E5)(75iV)7 + • • • , 


(B.9) 


where the ellipses stand for other Lorentz structures that do not contribute to In¬ 

serting everything into Eq. ([B.3|) we arrive at 


(Px) / dx 3 e-“ 3 ^-^)vf^'^^(a: 3 ) = P-x 


Vxxl 


dt e 


-iP-xx^t 


Vi 
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Vx xz j dt - V 2 ) 


j + V^3 + V"4 + 21 ^ 5 ) • (B.IO) 

In order to solve this equation we expand both sides at short distances and obtain the mo¬ 
ments of with respect to 0:3 expressed through moments of the distribution amplitudes 
dehned as = J Vx x^Vi{xi). One hnds 

1 


dx3a;3 Vi^^''^(a;3) = 


(n -M)(n -h 2 ) 


(-21/1 + 1/3 + 14 + 21/5) 


(B.ll) 


{n + 3 )v^WO+ 2 ) _ L(^Vi - - (Idi + 


{n -h l)(n 3) [ 

up to contributions of multiparton distribution amplitudes with extra gluons that have been 
neglected. 

The analysis of the n-quark contribution is performed in a similar way. We consider the 
matrix element 


a;“ [M*(0)C'7„n^(a;)] ^^(0) \P) = -x^ 


Vi 


Vf ) Pa (75iV), 


+ V2MP„ + V3M ( 7 a 75 iV)., + V,M^X^ (75iV)^ + 

and hnd repeating the same steps that lead to Eq. (E3): 

(0|e«‘[ti‘(0)CA'(i)]4(0)|F) = (0| [£«* [«‘(0)C;!;tl^(a:)]<i;(0)],JP) 


(B.12) 


™2 pi a2 

+T i di5;^{0k«‘7(0)CA'(to)]<(0)l^) 

= (0| [£«‘ [u‘(iS)Cfu’(x)] d‘(0)]„|P) + gluons, (B.13) 

the only difference being that the term corresponding to a total translation does not arise 
in this case. For the moments with respect to X 2 we get 


dx 2 a ;2 


{n -h l)(n -h 3) 
1 


(n + 3 )v^O) 0 + 2 ) - i(I/i - 1 / 2 )^“^^”+^^ 

- 0 , , 0,3 (-2V1 + E3 + 1/4 + 2E5)(“)("+2)1 _ ( 3 ^ 14 ) 

(n -M)(n 2) L J 

The corresponding expressions in the momentum fraction space are easily obtained by in¬ 
serting the conformal expansions for Vi,... ,Ve and inverting the moment equations. The 
result reads 


v"‘“’fe) = ||(A.c: + /«7), 

J2 


(B,16) 


= ||(A,7 + /«7) 
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with 


= -(1 - X2f [13 - 20/f + 3X2 + 10/r(l - 3x2)] , 
c; = (1 - X 2 )^ [113 + 495X2 - 552x^ + 10A“(-1 + 3 x 2 ) + 2Vi^^(113 - 951x2 + 828x^)] , 
Cf = -(1 - X 3 ) [11 + 131x3 - 169x2 + 63x^ - 30/f(3 + llxg - ITx^ + Tx^)] 
- 12 ( 3 - 10 / 1 ") ln[x3], 

Cf = -(1 - X 3 ) [1441 + 505x3 - 3371x^ + 3405x^ - 1104x^ 

-24i/i"(207 - 3 x 3 - 368x2 ^ 4^2x3 - 138x^)] - 12(73 - 2201/") ln[x 3 ]. (B.16) 
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